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Abstract 



Qh! We study the sensitivity of the HELLAZ and Borexino solar neu- 



trino experiments on discriminating the neutrino species i^e, i^e, ^fi,T, 
i^fi^Ty and using the difference in the recoil electron kinetic energy 
spectra in elastic neutrino-electron scattering. We find that one can 
^ ' observe a non-vanishing fo^,- component in the solar neutrino flux, 

' especially when the z^e survival probability is low. Also, if the data 

turn out to be consistent with z^e ^ '^^i,t oscillations, a component 
can be excluded effectively. 



1 Introduction 



The flux of solar neutrinos was first measured in the Homestake mine (see |jT| 
and references therein) over thirty years ago. Since then, it was reahzed that 
the measured flux was significantly suppressed with respect to theoretical 
predictions. More recently, a handful of different experiments have also 
succeeded in measuring the solar neutrino flux 0, |^, ^, |^. All experiments 
measure a neutrino flux which is significantly suppressed with respect to the 
theoretical predictions of the most recent version of the Standard Solar Model 
(SSM) 0. This thirty year old problem is what is referred to as the "solar 
neutrino puzzle." 

There are different types of solutions to the solar neutrino puzzle. At first 
sight, it appears natural to suspect that the SSM predictions for the solar neu- 
trino flux are slightly off, and/or that the experiments have underestimated 
their systematic effects, given that detailed models of the Sun and neutrino 
experiments are highly non-trivial. However, SSM independent analyses of 
the neutrino data (see for a particularly nice and simple example), together 
with independent experimental evidence in favour of the SSM Q, seem to 
indicate that the above solution to the puzzle is strongly disfavoured. 

The best solution to the solar neutrino puzzle involves extending the 
Standard Model of particle physics by assuming that the neutrinos have 
mass and that they mix, i.e., neutrino mass eigenstates are different from 
neutrino weak eigenstates. This possibility has become particularly natural 
in light of the recent strong evidence for i/^ oscillations from the atmospheric 
neutrino data at SuperKamiokande [||. 

Nonetheless, in order to firmly establish that the solution to the solar 
neutrino puzzle involves physics beyond the Standard Model, it is necessary 
to come up with SSM independent, robust experimental evidence for, e.g., 
solar neutrino oscillations. Indeed, these "Smoking Gun" signatures of solar 
i^e ^ t'other osciUatious are among the present goals of the SuperKamiokande 
experiment, via the measurement of the day-night asymmetry of the solar 



neutrino data M and the recoil electron energy spectrum |10], and the SNO 



experiment |jTI[ , via the measurement of the charged to neutral current ratio. 



the day-night asymmetry of the data, and the recoil electron kinetic energy 
spectrum. 

Other goals of this and the next generation of neutrino experiments are, 
if solar neutrino oscillations are established, to determine neutrino oscillation 
modes and measure masses and mixing angles. The current data allow for 
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different v^, oscillation modes and a handful of disconnected regions in the 
mass-mixing-angle parameter space (see [|12], |13| for two-flavour analyses and 
also [0 for an extension to the "dark side" of the parameter space). 

Experiments dedicated to measuring the flux of low-energy solar neutrinos 
[Ey = 0(100 — 1000) keV) are going to be extremely useful, and perhaps 
crucial, in order to fully solve the solar neutrino puzzle. It was recently shown 
that future experiments (Borexino and, perhaps, KamLAND [|l^]) dedi- 
cated to measuring the flux of ^Be neutrinos (produced by ^Be+e~ -^"^Li+Ue 
inside the Sun) should be able to establish or exclude the "just-so" solution 



I2| to the solar neutrino puzzle via the study of the seasonal variations of 



the neutrino flux , and establish or exclude the LOW MSW solution |T3 
via the study of the zenith angle dependence of their data |TB[. Furthermore, 
the measurement of a sizable ''Be neutrino flux would significantly disfavour 
the SMA MSW |TB| solution to the solar neutrino puzzle, especially in the 
case of Ue z/^ oscillations (where Ug is a sterile neutrino, i.e., a standard 
model singlet), and significantly constrain the SSM independent analysis, 
which require the flux of ^Be neutrinos to be virtually absent Finally, we 
have shown that, in the advent that the background rates at Borexino 
and/or KamLAND are exceptionally low, it should be possible to measure a 
nonzero component of i/^^r in the solar neutrino flux by analysing the recoil 
electron kinetic energy spectrum. 

Another exciting possibility is that of measuring the "fundamental" pp- 
neutrinos, which are produced in the interior of the Sun by proton-proton 
fusion {p + p ^^H+e+ + Ug) in a real time experiment. Future experiments, 
such as HELLAZ, HERON, LENS, etc (see for an overview) are being 
designed to do just that. The flux of pp- neutrinos is particularly constrained 
by the photon flux i.e., the Sun's luminosity, which, of course, is very 
well measured on the Earth. These lowest energy solar neutrinos {Ei, ^ 
420 keV) are not only the most abundant ones, but also have the best known 
flux. Their energy spectrum is also very well known, since it is dictated 
by the particularly well studied p + p nuclear fusion reaction. Among these 
proposed experiments, HELLAZ ||2T| will be able to determine the incoming 
neutrino energy in an event-by-event basis and have the unique opportunity 
of studying the solar neutrino spectrum and the recoil electron kinetic energy 
spectrum separately. Similar to what was shown for ''Be neutrinos [0, the 
authors of showed that HELLAZ may be able to measure a nonzero 
component of u^ ^- in the solar pp-neutrino flux by analysing the recoil electron 
kinetic energy spectrum independent of the SSM prediction for the solar 
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Table 1: Coefficients A, B in Eq. for different neutrino species. 



species 


A 


B 




sin^ Qw ^ \ 


sin^ Qy^ 




sin^ Q-w 


sin^ 6w + I 




sin^ Qw — \ 


sin^ 9w 




sin^ Qw 


sin^ 9w — \ 



neutrino ffux. 

In this paper we extend the analysis done in [19|, and study the ffavour 
composition of the flux of pp and ^Be neutrinos using the recoil electron 
kinetic energy spectrum. In particular we will address the capability of future 
low-energy solar-neutrino experiments to see evidence for z/^ coming from 
the Sun, and, in light of such evidence, exclude more "exotic" oscillation 
scenarios, such as Ve ^ t^s or ^ t'any oscillations. 

Our presentation is organised as follows: Sec. 2 describes the flavour de- 
pendent recoil kinetic energy distribution of events at Borexino and HELLAZ. 
Sec. 3 presents the technique for determining the presence of z/^^^ coming 
from the Sun, independent of the SSM prediction for the neutrino flux. We 
present simulations for both Borexino and HELLAZ and show how such a 
determination can be improved once we take the SSM prediction for the 
neutrino flux into account. Sec. 4 describes how the same procedure can be 
used to exclude the presence of antineutrinos or sterile neutrinos in the solar 
neutrino flux. In Sec. 5, we conclude. 



2 Recoil Electron Kinetic Energy Spectrum 

In this section, we discuss the differences in the recoil electron kinetic energy 
spectra among different neutrino species. Low-energy solar neutrinos are 
detected via 'V"+e~ ^"z/" +e~ elastic scattering in the experiments which 
will be considered here. By 'V" in the previous sentence, one actually means 
any of Ug, z/^, 1/^-, z/g, u^, or z/,-. Because z/^ and u^- are indistinguishable as far 
as the reaction above is concerned, we will refer to both as u^. 

The kinetic energy distribution of the recoil electrons, for a given incoming 
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neutrino energy is very well known and given by |^ 



dT ~ n 



(2.1) 



where rUe is the electron mass, T is the kinetic energy of the recoil electron, 
and Gp is the Fermi constant. The parameters A and B are given in Table |1|. 
The sign difference in the term 1/2 is a consequence of the presence (absence) 
of VT-boson exchange and the interchange of A and B between neutrino 
and anti-neutrino cases is a consequence of the "handedness" of the weak 
interactions. Eq. (|2.1|) is a tree-level expression, but higher order corrections 
are known to be very small |2^ , especially for the neutrino energies of interest. 



and will be neglected throughout. 

Borexino (under construction) is an ultra-pure liquid scintillator tank 
which detects the scintillating light produced by the recoil electron absorbed 
by the medium. For more details see |T^, |T^. It is sensitive to recoil 
electron kinetic energies greater than 250 keV, and is therefore sensitive to 
the (almost) monochromatic ^Be neutrinos with E,y = 862 keV. The expected 
resolution for the kinetic energy measurement varies from roughly 12%, for 
T = Tmin = 0.25 MeV, to 7% for T = T^^^ = 0.66 MeV H. They expect 53 
events/day in the SSM (BP95) together with 19 background events/day with 
the anticipated radiopurity of the scintillator of 10~^^g/g for U/Th, 10~^^g/g 
for ^'^K, and ^^C/^^C= 10~^^ and no radon diffusion. It is remarkable, 
however, that the Miinchen group of Borexino achieved a radiopurity for 
an organic liquid (Phenyl-ortho-xylylethane) better than 1.0 x 10"^''g/g P3| ; 
this is an upper bound on the contamination, limited by the sensitivity of 
the neutron activation measurement and hence the actual radiopurity may be 
even better. In this paper, we ignore the background to the ''Be solar neutrino 
signal at Borexino. This is probably an overoptimistic assumption, but could 
be realised in future upgrades given the above-mentioned achievement. 

HELLAZ (proposed) is a large time projection chamber (TPC) filled with 
roughly 2000 m^ of cool helium gas (~ 6 tons at 5 atmos, 77 K), which 
serves as the target for iz-e scattering. The recoil electron propagates in the 
gas medium before being absorbed, leaving a track of ionization electrons. 
These are then collected, yielding information about the kinetic energy and 
the flight direction of the recoil electron. HELLAZ is sensitive to recoil 
kinetic energies greater than ~ 50 keV, and can therefore "see" most of the 
pp-neutrino spectrum. Most importantly, since not only the recoil kinetic 
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energy of scattered electrons is measured but also their direction, it is possible 
to reconstruct the incoming neutrino energy, given that the position of the 
Sun in the sky is known, via the simple kinematic relation 

_ 2cos^g 

'"'(l+m,/E,)2-cos2^' ^ ' 

where 9 is the recoil electron scattering angle with respect to the incoming 
neutrino direction in the laboratory frame. Incidently, from Eq. (|2.2| ) it 
is very easy to compute the maximum value of the recoil electron kinetic 
energy, T^ax = T{e = 0) = EJ{1 + me/{2E^)). HELLAZ expects to 
measure the recoil electron kinetic energy with a resolution which varies 
roughly from 2% to 4% and the incoming neutrino energy with a resolution 
which varies between 5% and 12% PTl. They expect around 7 events/day 



from pp neutrinos in the SSM with negligible background. The major sources 
of background at HELLAZ are radioactive impurities from ^^^Th and ^^^U 
in the structure of the TPC However, because of the detector's total event 
reconstruction capabilities (including directional information), very efficient 



background rejection schemes are possible (see |^ and references therein for 
further information). 

The issue we would like to concentrate on is whether the shapes of the 
recoil electron kinetic energy distributions for different (anti)neutrino species 
are statistically different at Borexino and HELLAZ. With this in mind. Fig. ^ 
depicts the normalised distribution of events at HELLAZ (left) and Borexino 
(right). In the case of Borexino, the data is binned into ten kinetic energy 
bins, between 250 keV and 650 keV. In the case of HELLAZ, the data is 
binned into 4 x 21 bins in E^, x T. The bins have a width of 50 keV in the 
E^ direction and central values of 245, 295, 345 and 395 keV, while in the 
T direction they have a width of 10 keV in the range from 50 to 260 keV. 
The bin sizes have been chosen such that they are roughly the same as the 
resolution of both detectors. In order to integrate over the incoming neutrino 



*In addition to pp-neutrinos, HELLAZ is also sensitive to ''Be neutrinos, as well as the 
pep-neutrinos and the neutrinos coming from the CNO-cycle. ''Be neutrinos can be clearly 
separated from pp-neutrinos, while the number of pep and CNO-cycle neutrino generated 
events is expected to be less than 10% that of pp-neutrinos. Borexino is sensitive to, 
in addition to ''Be neutrinos (with E^, — 862 keV), a fraction of the pep and the CNO- 
cycle neutrinos, which produce approximately 10% as many events as ^Be neutrinos. We 
assume throughout, for simplicity, that only pp (''Be) neutrinos are detected at HELLAZ 
(Borexino). 
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energy at HELLAZ, the (normalised) BP98 pp-spectrum presented at was 
used. 




Figure 1: Normalised recoil electron kinetic energy distributions, for each of 
the 4 neutrino energy bins (see text) at HELLAZ (left) and for ^Be neutrinos 
(right). 

Many important features of the recoil electron kinetic distributions are 
worthwhile to point out. First of all, it is quite clear that the spectrum 
produced by z/g-e is much steeper than all the other ones.[[] Second, the v^. 
and generated spectra have opposite slopes when the neutrino energy is 
small enough, while their shapes start to look more and more similar as the 
neutrino energy increases. Finally, the spectra produced by i/^-e and v^-e 
scattering are extremely similar, especially at very low energies. 



All of these features can be readily understood from Eq. (2J.). First, it is 



convenient to write the expression for the normalised recoil electron kinetic 



^Some of these features were pointed out in 
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energy distributions, 



dy ^ 



^ B 2 me 



(2.3) 



where y = T/E^, and is a normalisation constant, such that / ^ from 
Tmin to Tmax equals unity. 

For Ue, ^ ~ 3, while ^ ~ 1/3. For the situation is reversed, while for 
both i/^ and ^^/x, ^ ~ ^ ~ —1-0 Note, curiously enough, that the reason for 
the similarity between the z/^ and z/^ cases is simply due to the accidental 
fact that sin^ 9w is close to 1/4. | This similarity is even more pronounced at 



very low energies, when the ^ term dominates over the ^ and ^ terms. 

Keeping in mind that 0.59 ^ ^ 2.3 in the neutrino energy range 
of interest and < y < (1 + me/2Ey)^^ < 1, one may write approximate 
expressions 

(2.4) 

(f)p.^3(l-y)2-^l/, (2.5) 
+ + (2-6) 

In the limit ^ <C 1, all three distributions are quite different (see, e.g., 



Fig. 1(B) in |T3)- The z/g case is roughly flat, the Ve case ranges from 3 at 
?/ = OtoOat?/ = l and the z/^, z/^ case ranges from 3/2 at y = to 3/4 at 
y = l. 

For ^ ^ 1, things are slightly more complicated, but still easy to 
understand. For example, the slope of the distributions for small values 
of y are, up to normalisation factors, — — ^6 + and + — 2^ for 
z/e, Ve and z/^, z/^ respectively. It is then easy to note that the z/g slope is 
significantly more negative than the other two, and that, in the case of z/^, z/^ 
the slope is actually positive if E,^ is small enough. This is indeed what one 
observes in Fig. ^ 

As the incoming neutrino energy increases, the distributions generated 
by z/g and v^, look more and more similar. One hint of this behaviour is 



■fin this case, the normahsation constant N is negative. 

§The fact that there is a sign difference between and ga for muon-type 

(anti)neutrinos is irrelevant, since these coefficients either appear as squares or as the 
product gLgR- 
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that the slope of the z/g case increases (decreases in absolute value), while 
the slope of the i/^, z/^ decreases. One can easily estimate that for 0.9 MeV< 
El, < 1.0 MeV the shapes of the Ue and z/^ induced recoil kinetic energy 
distributions are most similar. Indeed, for ^Be neutrino energies, one can 
already note that the difference between the z/g and the z/^ cases is similar to 
the difference between the z/^ and the cases. 

3 Measuring a z/^^ Component in the Solar 
Neutrino Flux 

In this section, we address the question whether the shapes of the recoil 
electron kinetic energy distributions presented in Sec. 2 are statistically 
different at Borexino or HELLAZ. In the affirmative case, there is hope that 
one may be sensitive to a "contamination" of other neutrino types in the solar 
neutrino flux by analysing the shape of the recoil kinetic energy spectrum. We 
consider this an "appearance experiment" of the "wrong" types of neutrinos 
from the Sun. In this section we will only consider the case of z/g ^ z/^ 
oscillations. 

In the advent of neutrino oscillations, a mixture of different neutrino 
weak eigenstates reaches the Earth. Given an electron-type neutrino survival 
probability Pee, a fraction Pee of all the neutrinos arriving at the detector 
are z/g, while a fraction 1 — Pee are z/^. The recoil electron kinetic energy 
distribution will, therefore, be given by 



da{T,E,] 
dT 




(3.7) 



Note that, in general. Pee is a function of the neutrino oscillation parame- 
ters (the mass-squared differences of the neutrino mass eigenstates and the 
neutrino mixing angles) and the neutrino energy. 

We simulate "data" at Borexino and HELLAZ for different values of 
Pee- We use the distributions presented in Sec. 2, while the flux of pp and 
^Be neutrinos are taken from the SSM In the case of Borexino, the 
energy dependence of Pee is irrelevant, given the monochromatic nature of 
^Be neutrinos. In the case of HELLAZ, we assume that Pee is constant inside 
each individual neutrino energy bin. Following the central idea presented in 
we perform a fit to the "data" using a linear combination of z/e-e 
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scattering and i/^-e scattering with arbitrary coefficients, 




(3.8) 



i.e., we perform a two parameter (Ce and C^) fit to the data. This mea- 
surement procedure is independent of the SSM prediction for the neutrino 
flux. Therefore, if a nonzero coefficient of the u^-e scattering distribution is 
measured, one can claim to have detected evidence for neutrinos other than 
Ue coming from the Sun. This "appearance" result certainly qualifies as a 
smoking gun signature for neutrino oscillations. 

Fig. 1^ (long, thin error bars) depicts the measured value of 1— Pee = c'^c 
in each of the neutrino energy bins defined in Sec. 2 as a function of the input 
value of Pee, for 5 years of simulated HELLAZ data. As was mentioned in 
the previous paragraph, the relevant information one should obtain from the 
plot is if the measured value of 1 — Pee oc is statistically different from 
zero. 

Fig. 1^ (right, long, thin error bars) depicts the measured value of 1 — 
as a function of the input value of Pee, for two years of Borexino running 
This is just a repetition of Fig. 2(A) in [|T^.[] Fig. ^ (left, long, thin error 



ee 



bars) depicts the result obtained at HELLAZ if all energy bins are used in 
the "data" analysis. This result is only meaningful if Pee is roughly constant 
for neutrino energies ranging from from 220 keV to 420 keV. This happens 
to be the case for most of the currently preferred regions of the two-neutrino 
oscillation parameter space, especially LMA, LOW and VAC solutions (see, 
e.g., [|12|)-C] Clearly, the significance of the measurement is better than the 



one obtained for individual energy bins (Fig. 

Next, the same analysis as above is repeated, except that the SSM predic- 
tion for the solar neutrino flux is included in the analysis. An uncertainty 
of 20% (5%) was assumed for the ^Be (pp) neutrino flux. The theoretical error 
was considered Gaussian for simplicity.^] Note that the uncertainties assumed 

*In a different variable, P = f — Pee, was used. Both results are, of course, 
equivalent. 

tPor the SMA solution, there is a sharp drop in Pee at E^, ~ 0.4 MeV, and the three 
lower bins can be combined without any problem. At HELLAZ, this will show up in the 
data, as the spectrum differs from the expected pp-neutrino spectrum, and hence is 
not a concern. 

■I- This procedure follows the one used in The readers are referred to this article 

for details. 
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Figure 2: The "measured" value of 1 — Pee as a function of the input value of 
Pee for each of the 4 neutrino energy bins (see text), after 5 years of HELLAZ 
running. The long, thin error bars correspond to model-independent analyses 
based only on the electron recoil energy spectrum shape, while the short, 
thick ones correspond to analyses which include the SSM prediction for the 
solar pp-neutrino flux, with an (inflated) uncertainty of 5%. 

here are inflated with respect to the ones quoted in the SSM calculations 
(9% and 1%, respectively), in order to render the results very conservative. 
Since the rates are very high both at Borexino and HELLAZ, the error 
bars are dominated by the uncertainties in the fluxes and hence they can 
be approximately scaled according to the assigned flux uncertainties. 

The results are presented in Figs. ^ and ^ (short, thick error bars). The 
significance of the measured value of 1 — Pee improves significantly, especially 
at HELLAZ, because of the small assigned uncertainty on the pp- neutrino 
flux. After five years of HELLAZ running, for example, one should be be 
able to determine a 1-sigma-away-from-zero z/^ component in the pp-neutrino 
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Figure 3: The "measured" value of 1 — Pge as a function of the input value of 
Pee, after 5 years of HELLAZ (left) and 2 years of Borexino running (right). 
Pee is assumed to be constant for E,^ =220-420 keV in the case of HELLAZ. 
The long, thin error bars correspond to model-independent analyses based 
only on the electron recoil energy spectrum shape, while the short, thick 
ones correspond to analyses which include the SSM prediction for the solar 
pp (jBe) neutrino flux, with an (inflated) uncertainty of 5% (20%). 



flux even for Pee ~ 0.9. It is also noteworthy that in the case of the SMA 
MSW solution to the solar neutrino puzzle Pee ~ for ^Be neutrinos, in 
which case a 4-sigma-away-from-zero evidence for z/^ in the solar neutrino flux 
can be established in only two years of Borexino running! It is important 
to emphasise that using SSM predictions for the solar neutrino flux is a 
reasonable thing to do, especially for pp-neutrinos. As mentioned before, the 
flux of pp-neutrinos is very well known because it is tightly related to the 
flux of light coming from the Sun. It is, therefore, the neutrino flux which is 
least sensitive to detailed modelling of the Sun's innards. 

Some comments are in order. First, only statistical uncertainties were 
considered, and there are no background events in our "data." As discussed in 
Sec. 0, the assumption of a negligible background rate seems less than realistic 
at Borexino, but may be possible in future upgrades. It may, however, be 
a fair assumption in the case of HELLAZ. If the real experimental data 
contains a sizable number of background events, it is necessary to either 
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subtract the background in a bin-by-bin basis or to somehow model the 
recoil kinetic energy distribution produced by background events. Analysing 
either of these procedures, however, is beyond the scope of this paper. 

Second, the analysis which does not include the SSM flux predictions 
is completely model-independent (the only assumption being the electron 
recoil spectrum as predicted by the standard electroweak theory), while the 
one which includes the SSM flux predictions is model-dependent. Obviously, 
one obtains a much better determination of Pee with the additional input of 
the SSM flux predictions. For establishing the "wrong neutrino component" 
in the solar neutrino flux as a smoking gun signature of the solar neutrino 
oscillations, the former approach is desired. However, for the purpose of 
determining the oscillation parameters, the energy dependence of the survival 
probability, and excluding other neutrino oscillation modes, such as Ve ^ 



s 



or z/g Ve,ii,T (as will be discussed in Sec. |), it is reasonable to include the 
SSM predictions in the analysis. 

Finally, we point out that the results we obtained for HELLAZ are 



similar to the ones obtained by J. Seguinot et al [22|. Indeed, we chose 



neutrino energy bins at HELLAZ which coincide with the ones used in 



22| . They also perform two different analyses of their simulated data, one 
which is independent of the SSM prediction for the solar neutrino flux, and 
one which assumes the SSM prediction for the flux. However, their data 
analysis procedure is somewhat different, and they do not take the theoretical 
uncertainty of the solar neutrino flux prediction into account. 



4 Testing for the z/g ^ i^s or ^e,/i,r Hypotheses 

Although it is most natural to assume that electron-type neutrinos oscillate 
into some linear combination of muon-type and tau-type neutrinos, there is a 
logical possibility that electron-type neutrinos might oscillate into standard 
model singlet sterile neutrinos |2^, or, perhaps, into antineutrinos of all 
flavours^ (see |^ and references therein). In this section, we will address the 



issue of excluding these solar neutrino oscillation modes if the data collected 
at Borexino and HELLAZ are consistent with z/g ^ z/^ oscillations. 

One can already address these "exotic" oscillation modes with the current 
experimental data. The flux of electron-type anti-neutrinos from the Sun 



*The original neutrino oscillation paper by Bruno Pontecorvo pOtl did, after all, consider 
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is particularly constrained by the SuperKamiokande and the LSD experi- 
ments the 95% CL SuperKamiokande upper bound on the flux of 
from the Sun with energies ^ 6.5 MeV is < 1.8 x 10^ cm~^ s~^, or 
^v^l^ssM — 0.035, where ^ssm the SSM prediction for the neutrino 
flux. KamLAND, being a dedicated detector for z/g, will improve this bound 
further down to 0.1% of the SSM flux above reactor anti-neutrino energies 
(Ey ^ 8 MeV) after one year of running [1^. There are, however, scenarios 
in which, for energies below the SuperKamiokande threshold, the v^, ^ Ve 
mixing is quite large ( ||2^ and references therein) . Such a possibility can only 
be addressed by low-energy solar neutrino experiments. 

Below, we discuss the exclusion of electron-type neutrino oscillations into 
sterile neutrinos or into one of the antineutrino types separately at Borexino 
and HELLAZ experiments. 



4.1 Z^e ^ 



The Vf. ^ Ug oscillation mode is allowed by the analysis of current solar 
neutrino data |jl3| , even though in the case of the MSW solutions to the solar 



neutrino puzzle, only the equivalent of the SMA MSW solution exists at the 
confidence level |T^. It is curious to note that, in the case of atmospheric 



neutrinos, the Us hypothesis is currently somewhat disfavoured 

In the case of z/g <^ Vs oscillations, one expects the recoil electron ki- 
netic energy spectrum to be exactly the same as the one generated by z/g-e 
scattering, since Vg do not interact with electrons. The only effect of the 
neutrino oscillations would be to suppress the expected number of events, 
i.e., the hypothesis of v^. oscillations is identical to assuming that the 

solar neutrino flux is, somehow, suppressed. Therefore, we attempt to fit 
the "data" simulated according to Eq. ( |3.7| ) in Sec. 3 (remember that the 
"data" is consistent with z/g ^ oscillations) to the trial function Eq. ( p.8|) , 
where the piece which corresponds to vanishes identically. This is a one 
parameter fit to Cg. Note that the only discrimination against Vs is the 
recoil energy spectrum, because the rate can be always fitted with the free 
parameter Cg. The inclusion of the SSM flux prediction does not help in 
excluding the Ve J^s oscillation because the free parameter Ce makes the 
predicted flux irrelevant. [] Note that ^Be neutrinos are predicted to have 

^One can "discover" Vg by observing a nearly vanishing rate for ^Be neutrinos. In the 
case of ^ v^^t oscillations, neutral-current scattering guarantees at least (when P^e = 0) 
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almost completely oscillated to Ug for the (only available) SMA solution: 

P,, = O.OO9l°;g^H0]- 

Fig. ^ depicts the value of obtained when one attempts to fit the "data" 

to the sterile neutrino hypothesis, for 5 years of HELLAZ running (left) and 2 

years of Borexino running (right). In the case of HELLAZ, we have assumed 

that Pee is constant over the entire pp-neutrino energy range. The value 

of is determined using the philosophy employed in [0, and should be 

compared to Nuns — 1 {Nuns is the number of "data" bins). After 5 years 

of HELLAZ running, one should be able to exclude sterile neutrinos coming 

form the Sun at more than 99.9% confidence level (CL) if all electron-type 

neutrino have turned into muon-type neutrinos {Pee = 0). After 2 years of 

Borexino running, the sterile neutrino hypothesis is only ruled out, at best, 

at the 89% CL.0 The explanation for this is the fact that the recoil electron 

kinetic energy spectra are very different when one compares the Ue-e and the 

z/^-e scattering cases at very low energies, i.e., pp- neutrinos, and similar at 

O(MeV) energies, i.e., ^Be neutrinos, as discussed in Sec. 2. The exclusion 

CL decreases with increasing Pee, and sterile neutrinos are excluded after 5 

years of HELLAZ running only at the 77% CL for Pee = 0.4. 

4.2 Ue ^ TJ, Model-independent Fit 

In the case of Ve ^ i^e,fj, oscillations, we perform a two parameter fit to the 
"data" simulated as in Sec. 3 to a linear combination of the Ue-e and i^e,^i-G 
scattering recoil kinetic energy distributions. Fig. ^ depicts the value of 
obtained when such a fit is performed, for 5 years of Borexino and HELLAZ 
running. The value of is to be compared to A^bins~2 to determine exclusion 
confidence levels. As advertised in Sec. 2, the z/^-e and z/^-e scattering 
cases produce almost identical recoil kinetic energy spectra, and are almost 
undistiguishable at HELLAZ. At Borexino, however, the difference between 
z/^-e and z/^-e scattering is similar to the difference between the z/^-e and z/g- 
e cases, as mentioned in Sec. 2 (see Fig. |l[), and some discrimination seems 
possible. Furthermore, upon close inspection, one should note that the shape 
of the distribution produced due to z/e-e scattering is more similar to the z/^-e 
case than the z/^-e. Therefore, any z/^ component in the trial function makes 

21% of the SSM rate. For this purpose, one should rely on the SSM flux prediction. 

■''The situation does improve, of course, it more events are collected at Borexino. After 
5 years of Borexino running, for example, one can exclude sterile neutrinos for Pee ^0.1 
at more than 95% CL. 
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Figure 4: Minimum values as a function of the input value of Pee, obtained 
when one tries to fit the "data" (which is consistent with Ue t'^ oscillations 
and Pee = Pee(input)) with a z/g + i/^ distribution (see text), for 5 years of 
HELLAZ (left) and 2 years of Borexino running (right). The dotted lines 
indicate the 95%, 99% and 99.9% exclusion confidence levels. 



the value of larger, i.e., the minimum of is obtained when the coefficient 
of the z/^ component is zero.0 This is exactly what happens in the case of 
Ue i^e oscillations, at both experiments. Any Ue component in the flux 
makes the agreement between the theoretical function and the "data" worse, 
and again the best value of is obtained when the coefficient of the z/g-e 
scattering distribution is zero. 

One can see from Fig. |^ that, after 5 years of HELLAZ data, Ue coming 
from the Sun can be ruled out at more than 95% CL if Pee ~ 0.2, while 

^ i^fi oscillations are not constrained at all, even for Pee = 0. After 5 
years of Borexino data, both Ue ^ and Ve ^ oscillations are ruled out 
at more than 95% CL if Pee ^0.1. 

Even if the Ue iy hypothesis cannot be ruled out at some reasonable 
CL, one may still be able to place upper limits on the flux of anti-neutrinos 
coming from the Sun. In the case of z/g ^ Ve oscillations, it is straight- 
forward to place upper bounds on the flux of electron-type antineutrinos at 

§We only allow nonnegative coefficients of the distribution functions in the fits, for 
obvious reasons. 
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Figure 5: Minimum values as a function of the input value of Pee, obtained 
when fitting the "data" with a Ue + iy distribution (see text). The fit does 
not include the SSM prediction for the solar neutrino flux and is for 5 years 
of HELLAZ (left) and Borexino running (right). The dotted lines indicate 
the 95%, 99% and 99.9% exclusion confidence levels. 



both HELLAZ and Borexino. The 95% CL upper bounds on the Ue flux are 
depicted in Fig. p. Of course, for Pee ~ 0.2 (0.1) at HELLAZ (Borexino) 
the upper bound on the flux is meaningless, since the hypothesis of Ue is 
already ruled out at more than 95% CL. Note that the upper bounds on 
the antineutrino fluxes are normalised by the SSM prediction for the pp- 
neutrino flux = 5-94 x 10^° cm"^s"^ for the HELLAZ result, and the 

SSM prediction for the ''Be neutrino flux $55^/ = 4.8 x 10^ cm~^s~^ for 
the Borexino result. For comparison, the 95% CL SuperKamiokande upper 
bound is 3.5% of the SSM flux, while KamLAND will improve it to 0.1% after 
one year of running. However, both of them are only for the ^B neutrinos. 
Both the HELLAZ and (especially) the Borexino limits obtained from 5 years 
of "data" are competitive with the SuperKamiokande limit for lower energy 
neutrinos. 

In the case of z/g ^ z/^ oscillations the situation is more ambiguous, 
especially at HELLAZ.^ Not only are the minimum values of very small, 

^The same is true at Borexino if one assumes the SSM prediction of the total neutrino 
flux, as will be described later. 
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Figure 6: Upper limit on the flux of electron- type antineutrinos after 5 
years of HELLAZ (left) and Borexino (right) running. The upper limits 
are normalised by Standard Solar Model (SSM) prediction for the pp (''Be) 
neutrino flux at HELLAZ (Borexino). 

but in some cases (especially for small values of Pee) a zero z/^ flux is ruled 
out at more than 95% CL. In such cases, it seems that the reasonable thing 
to do is to measure the antineutrino flux, not determine upper limits! The 
only exception to this is the case Pee = 1, when the data looks exactly like 
the SSM prediction, without neutrino oscillations. Indeed, one can not only 
set upper limits on the antineutrino fluxes, but should also set limits to the 
flux. Such limits are presented in Table |^. 

It is worthwhile to comment that the information contained in Figs. |^, 
H, and in Table | is also valid for the case of any unknown source of solar 

Table 2: Model-independent 95% CL upper limits on the flux of solar 
muon-type neutrinos and antineutrinos, when the data after 5 years of 
HELLAZ/Borexino running is consistent with SSM predictions. 



Experiment 






HELLAZ 


1.18 


1.44 


Borexino 


1.62 


3.77 
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antineutrinos of the electron and the muon-types, not only neutrino oscilla- 
tions. This is because our "data" was analysed assuming that the total flux of 
solar neutrinos is unknown. We emphasise that Pee is the survival probability 
of electron neutrinos assuming that they oscillate into active neutrinos, i.e., 
^ ^11. oscillations. 



4.3 z/e V, SSM-dependent Fit 

Next, the same analysis can be repeated assuming that the solar neutrino 
flux is known within theoretical errors. Again, the value of is computed^] 
and compared with A^bins — 2 + 1 (the —2 corresponds to the two coefficients 
that are varied during the minimisation procedure and the +1 corresponds 
to the solar neutrino flux constraint). Fig. |^ depicts the minimised values 
of obtained with 5 years of HELLAZ (left) and Borexino (right) "data." 
The theoretical uncertainty on the pp (''Be) neutrino flux was taken to be 
2% (20%); we inflated the theoretical errors by roughly a factor of two from 
those in BP98. 

A comparison between Figs. ^ and ^ reveals that the exclusion confidence 
levels increase, sometimes significantly. For example, after 5 years of HEL- 
LAZ one can exclude Ve ^ Ve oscillations for virtually all values of Pee at more 
than 99.9% CL. This is mostly because the z/e has a total cross section which 
is significantly larger than i/^, and the oscillation v^. <-h> v^. cannot account for 
the large suppression in the event rate in the "data" (due to ^ z/^). Even 
the elusive Ug <-> z/^ case can be excluded at Borexino at more than 95% CL 
for Pee ^ 0.2. Note that at HELLAZ the ability to discriminate between 
z/^ and is still quite limited. It is worthwhile to comment that, unlike in 
the case of model-independent fits in Sec. [4.2| , the minimum value of is in 
general obtained for a nonzero coefficient of the u-e scattering distribution. 
The reason for this is that, even though the shape of the z/-e scattering recoil 
electron kinetic energy distribution is "more wrong," the contribution to the 
overall cross section is smaller than the z/e-e scattering case, and therefore one 
obtains values of the solar neutrino fiux which are closer to the theoretical 
ones by having a finite u component, decreasing the value of x^- 

Again, one may set upper limits on the antineutrino fiux. As before, there 
is some ambiguity with regard to setting upper limits for the fiux, because 

"This procedure follows the one used in The readers are referred to this article 

for details. 
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Figure 7: Minimum values as a function of the input value of Pee, obtained 
when fitting the "data" with a z/g + z/ distribution (see text). The fit assumes 
the SSM prediction for the solar neutrino flux with an (inflated) uncertainty 
of 2% (20%) for pp (''Be) neutrinos, for 5 years of HELLAZ (left) and 
Borexino running (right). The dotted lines indicate the 95%, 99% and 99.9% 
exclusion confidence levels. 



for almost all values of Pee 7^ 1 at both experiments a zero flux is excluded 
at more than 95% CL. On the other hand, the z/g ^ z/g oscillation hypothesis 
is almost completely ruled out by HELLAZ and the upper limits obtained 
at Borexino are not much better than the ones depicted in Fig. ^ For this 
reason, the equivalent of Fig. |] in the case at hand is not presented. 

Table |^ contains the obtained upper limits on the (anti) neutrino fluxes 
when Pee = 1, i.e., when the data agrees with the predictions of the SSM. 
Unlike the case of a free total flux analysis, the results presented in Table ^ 
assume that the total neutrino flux of neutrinos to be detected at HELLAZ 
and Borexino is the one predicted by the SSM, i.e., there is no "room" for 
other, yet unknown, low-energy solar neutrino sources. For this reason, of 
course, the bounds obtained are (in some cases) much more stringent. 

Finally, as argued before, we emphasise that fixing the value of the solar 
neutrino flux to its SSM value is a reasonable thing to do, especially for 
pp-neutrinos. In these "exclusion analyses" such a procedure is even more 
natural, especially if one keeps in mind that a theoretical hypothesis, i.e., 
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Table 3: 95% CL Upper limits on the flux of solar muon-type neutrinos and 
antineutrinos when the data after 5 years of HELLAZ/Borexino running is 
consistent with SSM predictions, assuming that the total pp (^Be) neutrino 
flux is the one predicted by the SSM, with 2% (20%) (inflated) uncertainty. 



Experiment 


'^uJ'^SSM 






HELLAZ 


0.14 


0.14 


0.13 


Borexino 


0.66 


0.68 


0.058 



i^e t'/i oscillations plus the SSM computed values for the solar neutrino 
flux, has been "conflrmed experimentally." 

5 Conclusions 

In order to unambiguously solve the solar neutrino puzzle, and to estab- 
lish the oscillations of solar neutrinos (if they occur), clear "smoking gun" 
signatures are required. Such signatures include a large day-night effect, 
anomalous seasonal variations, or an obvious distortion of the neutrino energy 
spectrum. Another unambiguous signature is a discrepancy between the 
number of charged current and neutral current events at SNO, which can be 
viewed as an "appearance" experiment of f^^,-. However, SNO can look for 
this "appearance" signature only for ^B neutrinos with > 6.5 MeV and 
hence similar studies for lower energy neutrinos such as ^Be and pp neutrinos, 
which arc less sensitive to details of the solar model, are important. 

We have argued in this paper that a careful analysis of the recoil kinetic 
energy spectrum at Borexino and HELLAZ serves as another "smoking gun" 
signature, in the sense that one may be able to infer, independent of the 
SSM prediction for the solar neutrino flux, the existence of i/^^t- coming from 
the Sun. It is worthwhile to emphasise that this is different from distortions 
in the incoming neutrino energy spectrum. In our case we are describing an 
"appearance" experiment, while the analysis of the neutrino energy spectrum 
is a (energy dependent) "disappearance" experiment. 

It is important to point out that, in our simple simulations, no background 
events were included. While this is probably an oversimpliflcation in the case 
of Borexino, it may well be a good approximation for HELLAZ. Moreover 
future upgrades of Borexino may reduce the background further according 
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to recent encouraging progress One should keep in mind that, even if 
the background rates are significant, the procedure we described may still 
be useful if the background can be successfully dealt with (one should not 
underestimate the ability and creativity of experimental physicists!). 

We have also included in the analysis the SSM prediction of the fiux 
of solar neutrinos. While the results obtained in this manner are model- 
dependent (they are not "smoking gun" signatures of neutrino oscillations), 
we found them very useful. This is a reasonable thing to do especially for 
j9j9-neutrinos, whose fiux is constrained well by the solar luminosity. This 
additional input makes the measurement of the oscillation probability more 
precise. 

Finally, we have argued that, if the data collected at Borexino and HEL- 
LAZ is consistent with z/^_^ oscillations, one can try to exclude other 

neutrino oscillation modes (z/g ^ and ^ ^e,iJ,,T) using the same pro- 
cedure or, at least, to set upper limits on the fiux of solar antineutrinos. 
Again we considered the possibility of constraining the solar neutrino fiux 
to the SSM predicted value. The main result we obtained is that Vf. ^ v^, 
oscillations can, in general, be excluded, while the i^e ^ ^^i.t case is much 
more elusive. Nonetheless, Borexino should be able to exclude Vf. ^ Vf^^r 
oscillations if the SMA MSW solution to the solar neutrino puzzle happens 
to be the correct one. 
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